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The role of dynamical cosmological Casimir efltect to phantom (constant w) and oscillating universe 
is discussed. It is shown explicitly that its role is not essential near to Big Rip singularity. However, 
the account of Casimir fluid makes the scale factor approach to Rip time to be faster. Rip time 
itself maybe changed too. 
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I. INTRODUCTION. 

The Dark Energy problem is one of the main task 
for theoretical physics of last years. Since it was 
suggested its existence in 1998, a lot of candidates 
have been proposed (for a review see [H) in order 
to give a natural explanation of what we observe. 
One of these proposals deals with dark fluids, which 
possess an inhomogeneous Equation of State (EoS) 
that may depend on time (see Ref. In some 

works, the EoS is assumed to be periodic in time 
such that it could explain the so-called coincidence 
problem (see and references therein). On the 

other hand, as some analysis of the observations sug- 
gest (Ref. 0), the EoS parameter w for dark energy 
might have already (or nearly in the future) crossed 
the phantom barrier (w < —1), which could imply 
some kind of future singularity (for a classification 
of future singularities see [^). In the current work 
we study a constant and an oscillating dark energy 
fluid that crosses this barrier and drives the Universe 
evolution to a Big Rip singularity. We analyze the 
implications that a Casimir term could produce on 
the evolution, as it was pointed in Ref.flQj and 
this Casimir effect could play an important role on 
the dynamics of the Universe. Then, the effects on 
the cosmological evolution and specially close to Big 
Rip singularity are explored. 



II. THE REVIEW OF THE OSCILLATING 
DARK ENERGY UNIVERSE 

First of all, we review the oscillating universe sce- 
nario (see and refs. therein). 
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We consider the universe filled with the ideal fluid 
(dark energy) where the equation of state (EoS) pa- 
rameter w is time-dependent: 



P = w{t)p . 
The conservation law is: 

p + 3H{p + p) = 
and the FRW equation 
3 
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where H = ^, with a{t) the scale factor. The com- 
bination of (0) and dl]) gives 



p + KV3{l + w{t))p^^^ =0 , 
which can be integrated as 



{Jdt{l + w{t))) 
Using ([3]), the Hubble rate may be found 
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H = h(t) 



3 / dt{l + w{t)) 



(4) 
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When w is a constant (not equal to phantom divide 
value), the standard expression is recovered 



H = 



3{l + w){t~ts) ■ 



(7) 



If J dt{l + w{t)) =0, H diverges, which corresponds 
to the Big Rip type singularity, which occurs when 
t = ts- In case -u; > — 1, for the expanding universe 
where H > 0, the cosmological time t should be 
restricted to be t > is- 

The situation where w{t) is time-dependent, pe- 
riodic function may be of physical interest, as it is 
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pointed in several works (see Q and @) . An inter- 
esting example is given by 



-1 + wq cos ujt . 



Here, wo,uj > 0. 
Then Eq.® gives 



H = 



2uj 

3 {wi + Wq ainujt) 



(8) 



(9) 



Here wi is a constant of the integration. When 
< Wq, the denominator can vanish, which corre- 
sponds to the Big Rip singularity. When \wi\ > wq, 
there does not occur such a singularity. Since 



H 



2aj^wo cos ujt 
3 {wi + Wq sinuit)'^ 



(10) 



when Wq cos ujt < {wq cos ujt > 0) , the universe 
lives in phantom (non-phantom) phase where H > 

(f<o). 



III. THE ACCOUNT OF CASIMIR FLUID 
IN OSCILLATING AND PHANTOM 
COSMOLOGY 

A simple and natural way of dealing with the 
Casimir effect in cosmology is to relate it to the sin- 
gle length parameter in the (fc = 0) theory, namely 
the scale factor a. It means effectively that we 
should put the Casimir energy Ec inversely propor- 
tional to a. This is in accordance with the basic 
property of the Casimir energy, viz. that it is a mea- 
sure of the stress in the region interior to the " shell" 
as compared with the unstressed region on the out- 
side. The effect is evidently largest in the beginning 
of the universe's evolution, when a is small. At late 
times, when a — > oo, the Casimir influence should be 
expected to fade away. As we have chosen a nondi- 
mensional, we shall introduce an auxiliary length L 
in the formalism. Thus we adopt in model in which 



C 

21^' 



(11) 



where C is a nondimensional constant. This is the 
same form as encountered for the case of a perfectly 
conducting shell. In the last-mentioned case, C was 
found to have the value 



C = 0.09235. 



(12) 



In the following we shall for definiteness assume C 
to be positive, corresponding to a repulsive Casimir 



force, though C will not necessarily be required to 
have the value p2)) . The repulsiveness is a char- 
acteristic feature of conducting shell Casimir the- 
ory, following from electromagnetic field theory un- 
der the assumption that dispersive short-range ef- 
fects are left out. Another assumption that we shall 
make, is that C is small compared with unity. This 
is physically reasonable, in view of the conventional 
feebleness of the Casimir force. The expression ((TT|) 
corresponds to a Casimir pressure 



-1 



dEr 
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Pc = 



ATi{Laf d{La) AnL'^a'^' 



(13) 



and leads consequently to a Casimir energy density 
Pc oc 

Repeating steps of [1(|] we have: 



Pc 



C 



Than, FRW equation of motion looks as 



3C 



(14) 



(15) 



as Eq.Q. Note that combination of two non- 
coupled fluids maybe considered as single inhomo- 
geneous fluid (see and @). 

The conservation law Eq.(l2]) is the same, because 
fluids do not couple: 



p + 3H{p + p) = 0. (16) 
Substituting p3|) into one gets: 



p + KV^{l + w{t))p^ip+^)^0. (17) 

In general, this equation is very difficult to be 
solved analytically. Although, by deriving the Fried- 
mann equation (jisp respect t and combining with 
the energy conservation equation (I16|) . it yields 



' + -ii + w{t))e^ 



-il~3w{t)) 



3GC 
2L4a4 



(18) 



where 9{t) — 3H{t). The case of viscous dark fluid 
[l3| is very similar. This equation describes the dy- 
namics of the Universe evolution for a given function 
w{t), and it is solved as a = a{t). It is easy to see 
that when the equation of state parameter w = 1/3, 
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the contribution of the Casimir effect to the dynam- 
ics of the Universe vanishes. Let us analyze first of 
all, a simple example when w(t) — wq is a constant. 
For this case, the equation (IT51) is rewritten as 



1 + 3wo 



-(l-3w;o) 



GC 



(19) 



2 a ' "^"""2X403 ■ 
And the general solution for this equation results: 

ada 
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where k is an integration constant. When we have 
a phantom field with wq << —1, the solution takes 
the form: 



a{t) = B 



1 + 



1 



l/l-3i«o 



here B ^ ^ 



gv kB(ts-t) _ I 

(21) 

Note that for t ^ tg, the scale factor 
given by ([2T|) diverges, such that the so-called Big 
Rip singularity takes place. On the other hand, if 
we neglect the Casimir contribution < C << 1, 
the solution for the scale factor yields 



a{t) ^ {ts - t)-s 



(22) 



which is the solution in absence of the Casimir term 
([T]). This solution grows slower than in the above 
case when the Casimir contribution is taken into ac- 
count (l?T|) . Also note that the Rip time tg for each 
case is different if the Casimir contribution is ne- 
glected in the second case, due to ts depends on the 
content from the Universe as it can be seen by the 
first Friedmann equation ts — tn= f°° -n — tv^^tttty, 

where to denotes the current time and r2'(a) is the 
cosmological parameter for the component i of the 
Universe 

Let us now to analyze the equation (|18p when w = 
w{t) is a periodic function dS]). By inserting ([5]) 
into (|18p . the equation turns much more compli- 
cated than the constant case studied above. As the 
Casimir contribution is assumed to be very small, we 
can deal the equation (jl8p by perturbation methods. 
Hence, the solution might be written as 



9{t) =00 + C0i + 0{C^) . 



(23) 



And by inserting ([23]) into (fT8|) . the zero and first 
order in perturbations can be separated, and it yields 
the following two equations: 

^0 + cosujtOo ~ , 



SG 1 

9i + Wo cos (jjtOoOi = — — — j- (4wo cos ut) 
2L^ 



(24) 



where ao(0 = exp (/ g^^^^^f^^;^^ 
lution for the zero order calculated in the previ- 
ous section, see ([9]). We suppose for simplicity 
wi Wo, such that a Big Rip singularity appears 
when 1 -|- sinujts = 0. As we are interested in the 
possible effects close to the singularity, we can ex- 
pand the trigonometric functions as series around ts . 
Then, the second equation in (|24|) takes the form: 



is the so- 



1 



3w2 ts — t 



3G , , 3wo , , , . . 16 

-(4H -{ts-t\)G 3^(t,-t) 



(25) 



whose solution is given by 

01 (i) = k{ts - t)^ + {ts - t)^g{llts - t) , 



where 



3GC 
2IJb 



4(4/1 + B) 
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Here A 



and B 



ie_ 

3u 



(26) 

Then, we can analyze 
the behavior from the contribution to the solution 
at first order ([26|) when the Universe is close to the 
singularity. When t ^ ts {x ^ oo), the function 
([^5]) goes to zero, what means that the Casimir con- 
tribution has no effect at the Rip time. 

Thus, we have demonstrated, that dynamical 
Casimir effect gives no essential contribution to 
(phantom oscillating) dark energy dynamics near to 
Rip singularity. 



IV. DISCUSSION 

We discussed the role of dynamical cosmological 
Casimir effect to phantom (constant w) and oscillat- 
ing universes. It is shown explicitly that its role is 
not essential near to Big Rip singularity. Moreover, 
the Rip time tg may change due to the account of 
Casimir fluid contribution. For oscillating universe 
approaching to finite-time singularity, the account 
of dynamical Casimir effect is made approximately, 
using perturbations technique. It is possible that 
only quantum gravity account fl3| may prevent the 
occurrence of Big Rip singularity. 
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